A detailed study on the range fluctuations of 10 n eV to 10 18 eV muons in 
water and the fluctuations of the Cherenkov light yields due to the 
accompanied cascade showers initiated by muons 
Parti: High Energy Muons 
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Abstract 

The validity of our Monte Carlo simulation procedure {the time sequential procedure) had been verified 
by the corresponding analytical procedure which is independent of our method methodologically Also, the 
results obtained by our procedure are compared with those obtained by the different Monte Carlo simulation 
procedure {the v cu t procedure) which has been exclusively utilized by different authors and the agreement 
between ours and theirs are found to be good. By utilizing our Monte Carlo procedure, the validity of which 
is guaranteed in two different procedures, we investigate not only the fluctuation of high energy muons 
themselves but also the fluctuation of the various quantities related to the energy losses by the muons, 
which arc difficult to obtain by the v cut procedure. Namely, we obtain the fluctuation on energy losses of 
muons in the present paper(Partf ). In a subsequent paper (Part2), based on the idea developed in Parti, 
we will examine the correlation between the primary energies of muons and their energy deposits and the 
correlation between the energy deposits of the muons concerned and their Cherenkov light yields in the KM3 
detector in detail. 

Keywords: High energy muon, Muon range fluctuation, Muon energy loss, Muon propagation simulation 



1. Introduction 

In the KM3 detectors deployed in the Antarctic, 
the ocean and lakes [l|, Q, H,1H i it is most important 
to decide the energies of high energy muons due 
to neutrino interactions from the experimental and 
technical point of view. In the energy determina- 
tion of high energy muons, there are two indispens- 
able problems to be carefully examined: the first is 
to examine the fluctuation of the traversed lengths 
of muons and the second is the fluctuation on the 
Cherenkov light yields initiated by both the muons 
concerned and their daughters' electromagnetic cas- 
cade showers. In the present paper (Parti), we 
studied the range fluctuation of high energy muons 
in detail and in a subsequent paper (Part2), we will 
study the fluctuation on the Cherenkov light yields 
of the electromagnetic cascade showers which influ- 
ence decisively over the energy determination of the 



high energy muons in detail. 

The studies on the range fluctuations of high en- 
ergy muons have a long history (H|-(l6j| even before 
the appearance of the high performance of elec- 
tronic computers, but the recent big progress of 
the computers makes it possible to treat fluctuation 
problems on high energy muons in a more quanti- 
tative manner using the Monte Carlo method. As 
far as the treatment of the range fluctuation of 
high energy muons by the Monte Carlo method is 
concerned, there exist two independent approaches. 
The first one is the procedure of splitting the radia- 
tive energy loss into two terms, a continuous "soft" 
term for v < v cut and a stochastic "hard" term for 
Vcut < v < 1 in the equation for energy loss (see, 



Eq.©) [17l-|23|. Hereafter, we will simply refer to 



it as the v cut procedure. The introduction of v cu t is 
for the purpose of applying the Monte Carlo meth- 
ods to the "hard" term and its value is carefully 
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chosen to save computational time. The other is the 
time sequential procedure developed by us in which 
the interaction points of the muons and their dissi- 
pated energies are directly determined 0, and 
here, all the processes are treated in the stochas- 
tic manner without the introduction of v cu t . These 
two methods are independent of each other, but are 
logically equivalent, giving the same results as for 
the muons' behaviors (see Figures |4] to 6). Here, it 
should be noticed that in KM3 detectors, the energy 
determination of high energy muons are carried out 
by measuring Chcrcnkov light yields which are pro- 
duced by not only the muons themselves but also 
by the electromagnetic cascade showers induced by 
the muons concerned. Then it should be empha- 
sized that the largest Cherenkov light yields are oc- 
cupied by those generated by electromagnetic cas- 
cade showers and the Cherenkov light yields pro- 
duced by the original muons can be neglected for 
muon energies above ~ 10 14 eV (See detail, in 
the subsequent paper [2 7]). These cascade showers 
are generated from the stochastic processes, such 
as bremsstrahlung, direct electron pair production 
and photo nuclear interaction which are initiated by 
muons concerned. In the present paper, we discuss 
the behaviors of high energy muons in watcr(ice) 
and the problems around Cherenkov light yields are 
treated in a subsequent paper. 

2. Range fluctuation of the (ultra-) high en- 
ergy muons and individual behavior of the 
muons 

2.1. The physical meaning of "no fluctuation" 

The average energy loss by high energy muons 
is usually described as, 



^ = a(E) + b(E)-E, 
ax 



(1) 



where a is the term due to ionization which is free 
from fluctuation and b is the term due to stochastic 
processes which may be origins of fluctuations. The 
latter is divided into three parts. Namely, 



b (E) = b b (E) + b d (E) + b n (E) , 



(2) 



where 6b, bd and b n are the corresponding terms due 
to bremsstrahlung, direct electron pair production 
and photo nuclear interaction, respectively. In the 
treatment of the average energy loss, each 6-term is 
defined as, 



do~i (v, E) 
dv 



minimum fractional energies due to their kincmat- 
ical limits. Here, N and A denote the Avogadro 
number and the Atomic number of the medium. 
The physical meaning of Eq.([T]) is that the muons 
concerned dissipate energy in the form of the aver- 
age uniquely, being defined by Eq.Q. 

Namely, fluctuations are not included in Eq.Q. 
In Figure [TJ we give the 6-terms due to different 
processes in water. In Figure [21 we give the rela- 
tion between a(E) and b(E). As the 6-terms are 
essentially of a stochastic character, it is seen from 
the figure that the stochastic processes become ef- 
fective above ~ lTeV(10 12 cV). Therefore, we must 
treat the muon's behavior in a stochastic manner 
above ~ 1 TeV. Below ~ ITeV we may treat muon's 
behavior in the non-stochastic manner. Then, the 
range of the muon is uniquely determined by Eq. (QJ . 



R 



E 



dE 

dE 

dx 



(4) 



where E m i n is the minimum energy for observa- 
tion and Eq is the primary energy of the muon. 
Throughout the present paper, E m i n denotes the 
minimum energy among the energies for observa- 
tion (Eobs) and is taken as lGeV (10 9 eV). Thus, 
R defined by Eq.Q gives the effective range of the 
muon without fluctuation. 

In Figure [3 we give the effective range defined 
by Eq.@ together with our average ranges of the 
muons in which the fluctuation effects are exactly 
taken into account. (see discussion in the later sec- 
tions). It should be noticed from the figure that 
the effective ranges without fluctuation arc differ- 
ent from the average range of the muons in which 
the fluctuation is considered, on which Lipari and 
Stanevl7| already pointed outQ. Really, the real 
average ranges are smaller than those of effective 
range beyond one standard deviation above 10 13 
eV as shown in Figure [3] 

2.2. Physical quantities with fluctuation 



In the v cut procedure, many authors |17|- [23[ di- 
vide all stochastic processes into two parts, namely, 
the continuous part and radiative part in their 
Monte Carlo simulation in order to consider the 
fluctuation in both the range and energies of the 



• dv. 



(3) 



1 What Lipari and Stanev pointed out had been confirmed 
by J.Nishimura in the analysis of intensity-depth relation of 
high energy muons. See page 111 Table 22. |25|| . 
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Figure 1: 6-terms due to different stochastic processes, a, 
b and c denote direct electron pair, bremsstrahlung and 
photo nuclear interaction, respectively. The lines denote 
ours, while symbols are due to Klimushin et.al. 
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Figure 2: The relation between a-term and b-term in water. 
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Figure 3: The range energy relation between the case with- 
out fluctuation and the case with fluctuation in standard 
rock. The uncertainty bar denote the standard deviation. 



muons and introduce v cut to save time for compu- 
tation, while we treat all stochastic processes as ex- 
actly as possible without the introduction of v cu t 



in the time sequential procedure [24j [26| and the 
essential structure for calculation is described in 
the Appendix. The validity of our Monte Carlo 
method had been already checked by the analytical 
method based on the diffusion equation [l2| which 
is methodologically independent of the Monte Carlo 
procedure. Furthermore, we check the validity of 
our method, comparing our results with the corre- 
sponding results by the v cut procedure. 

2.2.1. The comparison of our results with others 

The survival probability for high energy muon in 
our method is defined as, 



P [EohsiX, E p ) = 



^through 



{E P ) 



(5) 



where E p , A, and E b s , denote the energy of the 
primary muon, the point of observation and the 
minimum energy of the muon at the point of obser- 



vation, respectively. N s , 



(E p ) denotes the total 



sampling number of muons and N t hrough (E b s ,X) 
denotes the number of muons concerned with ener- 
gies above E b s which pass through the observation 
point X. 

We compare our results by the time sequential 
procedure with the different authors' results by 

the v cu t procedure in the following. Lipari and 

Stanev[17| give the survival probabilities as the 
functions of the depths for 1 TeV to 10 6 TeV in- 
cident muons in water and partly standard rock, 
the minimum energy of which is taken as 1 GeV. 
We obtain the corresponding results by the time 
sequential procedure and compare our results with 
Lipari and Stanev's in Figure 2J 

Also, Klimushin et al[20[ give the survival prob- 
abilities for primary energy of 10 13 eV to 3 x 10 16 
eV. Wc obtain the corresponding results to them 
and compare our results with their results in Fig- 
ure [5] 

Furthermore, Kudrvavtsev [22$ gives the energy 
spectrum of the muon due to primary energy of 2 
TeV at 3 km in water. We compare our results with 
his in Figure [51 

Thus, rather nice agreements between ours and 
the results by different authors shown in Figures 
S] to H2 are found to guarantee the validity of 
our method, taking into account the slight differ- 
ence in their differential cross sections such as the 
bremsstrahlung, direct electron pair production and 
photo nuclear interaction by different authors and 
the slight difference in their procedures, in addition 
to the confirmation of the validity of our method by 
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Figure 4: The comparison of our result with that of Lipari 
and Stanev fT^I . The survival probabilities of muons of en- 
ergy from 1 TeV to 10 6 TeV. The numerical figures attached 
each curve denote the primary energies. Curves labels corre- 
spond to following set of primay energies of muon: (l)lTeV, 
(2)10TeV, (3)10 2 TeV, (4)10 3 TeV, (5)10 4 TeV, (6)10 5 TeV, 
(7)10 6 TeV. Symbols are due to Lipari and Stanev and curves 
arc due to ours. 




Figure 5: The comparison of our results with that of 
Klimushin et al l20l . The continuous lines are obtained by 
us, while symbols are readout from those by Klimushin et 
al for primary energies from 10 13 eV to 3 X 10 16 eV. The 
numerical figures attached each curve denote the threshold 
energy is 10 TeV. Curves labels correspond to following set 
of depths: (1)1. 15km, (2)3.45km, (3)8. 05km, (4)12. 65km, 
(5)17.25km, (6)21. 39km. 



the analytical method. Here, we briefly mention to 
the essential difference between the v cu t procedure 
and the time sequential procedure. In the v cu t pro- 
cedure, one writes down the fundamental equation 
in the following, 



dE 
dx 



dE 
dx 



+ 



soft 



dE 



dx 



Figure 6: The comparison of our results with that of 
Kudrvavtsev [22Tl . Energy spectrum of muons with initial en- 
ergy of 2 TeV at 3 km. The line denotes ours, while symbols 
are due to Kudryavtsev 
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1 da(v,E) 

dv ■ v , (6) 



hard 



As Eq.© expresses the average energy loss of the 
muons, this method gives the same result concern- 
ing the energy loss of the muons concerned, irre- 
spective of the v CU f . In Eq.©, the Monte Carlo 
procedure is carried out in the second term and 
the difference in v cu t gives the different shape of 
the energy losses. The only and essential differ- 
ence between the v cu t procedure and the time se- 
quential procedure lies in that these two methods 
may give different shapes of the emitted energies of 
the muons even if their total dissipated energies are 
same. The difference in the shapes of the emitted 
energies between these two methods greatly influ- 
ences their shapes of the Cherenkov light yields and 
their magnitude, which, in turn, may give the dif- 
ferent estimation on the muon's energies concerned. 
We will examine carefully Cherenkov light yields 
problems in detail in a subsequent papcr(Part2) 

2.2.2. Survival probabilities and their v cut energy 
spectra at different observable depths 
In Figures [7] to [5J we give the survival probabil- 
ities for different observable energies with primary 
energies of 10 12 eV, 10 15 eV and 10 18 eV, respectively. 
In Figure[?l we give the observation energies 10 9 eV, 
10 10 eV and 10 11 eV, respectively. In Figure [5J we 
give them 10 9 eV, 10 10 eV, 10 n eV, 10 12 eV, 10 13 eV 
and 10 14 eV, respectively. In Figure[9j we give them, 
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Figure 7: The survival probabilities for 10 12 eV muon. 
Curves labels correspod to following set of cutoff energies: 
(l)10 9 eV, (2)10 10 eV, (3)10 11 eV. 
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Figure 8: The survival probabilities for 10 eV muon. 
Curves labels correspod to following set of cutoff energies: 
from (l)10 9 cV to (6)10 14 eV. 
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Figure 9: The survival probabilities for 10 18 eV muon. 
Curves labels correspod to following set of cutoff energies: 
from (l)10 9 cV to (9)10 17 eV. 



10 9 eV to 10 17 eV, respectively. Each sampling num- 
bers in Figure [7] to [9] are 100,000. It is seen from 



Figure 10: Energy spectrum in water at the different depths, 
initiated by 10 12 eV muons. 
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Figure 11: Energy spectrum in water at the different depths, 
initiated by 10 15 eV muons. 
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Figure 12: Energy spectrum in water at the different depths, 
initiated by 10 18 eV muons. 



the figures that the survival probabilities do not 
decrease remarkably as their primary energies in- 
crease. 

In Figures [10] to [121 we give the differential en- 
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Figure 13: Range distributions for 10 12 eV to 10 18 eV muons 
in water. The minimum observation energies are taken as 
10 9 eV. Each sampling number is 100,000. 



ergy spectra for primary energies, 10 12 eV, 10 15 eV 
and 10 18 eV at different depths, respectively. It 
is seen from the figures that the energy spectrum 
at the initial stage are of delta-function and they 
shift as the mountain-like deforming their shape 
in the intermediate stage and finally, they disap- 
pear in deeper depth for their threshold energy (1 
GeV). Each sampling number in Figures 10 to 12 is 
100,000. 

2.2.3. Range Distributions of Muons 

All processes, such as bremsstrahlung, direct 
electron pair production and photo nuclear inter- 
action are stochastic and, therefore, one cannot 
neglect their fluctuation essentially. Muons pene- 
trate through matter as the results of the competi- 
tion effect between bremsstrahlung, direct electron 
pair production and photo nuclear interaction. The 
logic on their penetration is described in the Ap- 
pendix. 

In Figure [T31 we give P (R; E ), the probabilities 
for the range distribution in water with primary en- 
ergies, 10 13 eV, 10 15 eV and 10 18 eV in water whose 
minimum energy is 10 9 eV, respectively. It is clear 
from the figures that the width of the range dis- 
tribution increases rapidly, as their primary energy 
increases. Also, as the primary energy decreases, 
the width of range distribution becomes narrower 
and approaches to a delta function-type, the limit of 
which denotes no fluctuation. It is interesting that 
the range distribution can be well approximated 
as the normal distribution above ~ 10 14 eV where 
the total Cherenkov light yields comes almost from 
the accompanied electromagnetic cascade showers 
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Figure 14: Hypothetical range distributions in water for 
10 13 eV muons together with the real range distribution. 
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Figure 15: Hypothetical range distributions in water for 
10 eV muons together with the real range distribution. 
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Figure 16: Hypothetical range distributions in water for 
10 18 eV muons together with the real range distribution. 



due to direct electron pair, bremsstrahlung, and 
photo nuclear interaction initiated by the muon(see, 



G 



Table 1: The average values, the standard deviations and the 
relative variances of the range distributions of muons from 
lO^eV to 10 18 cV in water. 
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Part2) and they are given as, 

where Eq, < R > and a are primary energy, the av- 
erage value of ranges and the standard deviations, 
respectively. Their average ranges, standard devi- 
ations and relative variances (standard deviations 
divided by averages) in water are given in Table [TJ 
Also, it is interesting that their relative variances 
decrease slightly as their primary energies increase. 

In order to examine each characteristic of 
stochastic process, such as the bremsstrahlung, di- 
rect electron pair production and photo nuclear in- 
teraction, we compare the hypothetical range distri- 
bution in which a specified stochastic process only 
is assumed to occur with the real distribution in 
which all stochastic processes occur as the result of 
the competition effect among these processes. We 
compare the real range distributions with the hy- 
pothetical range distribution in Figures Q3] to [T5] 

Here, the symbol d in these figures means the 
hypothetical range distribution in which only di- 
rect electron pair production is taken into account 
and the bremsstrahlung and nuclear interaction are 
neglected. The symbols b and n have similar mean- 
ing to that of d. The symbol t means the range 
distribution in which all interactions are taken into 
account (The turc distribution). From the shapes 
of the distributions and their maximum frequencies 
for different stochastic processes in Figures Q3] to 
ll~6l it is clear that energy losses in the direct elec- 
tron pair production are of small fluctuation, while 
both the bremsstrahlung and photo nuclear interac- 
tion are of bigger fluctuation and the fluctuation in 
photo nuclear interaction becomes remarkable when 
compared with bremsstrahlung as primary energy 



increases. 

2.2.4. The diversity of individual muon behavior 

In Figures [T7] to [23 we show the diversities of 
muon behaviors for the same primary muon ener- 
gies with regard to their ranges (or their energy 
losses). Also, in Table [2j we summarize the charac- 
teristics of the events in Figure [T7] to [37] In these 
figures, all interaction points due to the processes of 
bremsstrahlung, direct electron pair production and 
photo nuclear interaction and all energy losses at re- 
spective points due to these processes are recorded. 

In order to clarify the diversities among the real 
range distributions (or real energy loss distribu- 
tions), we show the muons with the shortest range, 
the muons with longest range and the muons with 
average like range. The average-like range is defined 
as the nearest range to the average range which is 
calculated from the total sampled ranges. 

In Figures [T7J to [T9j we give the characteris- 
tic behaviors with the shortest range, the average- 
like range and the longest range as a function of 
the energy loss and the depths traversed for the 
same energy of 10 12 cV in watcid. In the Fig- 
ures [TjJ to [T9l we utilize the same scale in depth to 
clarify the diverse behaviors by the same incident 
energies, namely, those with the shortest range, 
with the average like-range and with the longest 
range, respectively. In figures, the abscissa denotes 
the depths where the interactions concerned occur. 
The 'needles' with different colors denote the en- 
ergy losses due to direct electron pair (green, d), 
bremsstrahlung (red, b) and photo nuclear interac- 
tion (blue, n), respectively The ordinate shows the 
ratio of the energy loss concerned to the incident 
energy Namely, 'needles' denote the energy losses 
due to the interaction concerned and line graphs 
(black curve) denote the ratios of the muon ener- 
gies at the depths concerned to their incident ener- 
gies. The abrupt changes in them denote the catas- 
trophic energy losses for muons (see Figures [20] and 

It is seen from figures and Table [5] that there 
is not so big difference between the case with the 
shortest range and one with the longest range. In 
the case with the shortest range fFigurc [T7|) . we find 
two catastrophic energy losses (at ^910 meters and 



2 As for the characteristic behaviors of high energy muons 
which arc shown Figures 1171 to 1191 we suggest the readers 
that they look at the pictures with colors in the WEB page. 
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Figure 17: The energy losses with the shortest range for 10 12 eV muons. 
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Figure 18: The energy losses with the average-like range for 10 12 eV muons. 
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Figure 19: The energy losses with the longest range for 10 12 eV muons. (For interpretation of the references to color in legend 
of these Figures 17 to 20, the reader is referred to the web version of this article.) 



~1870 meters) due to brcmsstrahlung play the im- 
portant role in the range. In the case with average- 
like range (Figure ITS)) , we can find one catastrophic 
energy loss due to bremsstrahlung at ^1.48 kilo- 
meter. However, in the case of the longest range 
(Figure [IT)]) , we cannot find the catastrophic en- 
ergy losses due to bremsstrahlung and, instead, we 
can find that almost energy losses are due to many 
number(^300) of direct electron pair production 
events with smaller energies. It is seen from Figure 
[T3] that the fluctuation effect in the range distribu- 
tion is not so effective for 10 12 eV muon and their 
different behaviors with the different ranges are not 
so impressive. 

Figure [20] shows the same in Figure [21] in differ- 



ent scale. The shortest range, ^940 meter (Figure 
[2U]) . is far shorter compared with the longest one, 
^35.0 kilometers fFigure l23"|) . It is seen from Figure 
l20l and the Tables that brcmsstrahlung plays a de- 
cisive role as the cause of catastrophic energy loss, 
(-96.5% of total energy at ~450 meters). 86.6% 
of the total energy is lost by 2 bremsstrahlungs, 
13.4% by 367 direct electron pair productions and 
5.9 x 10~ 3 % by 1 photo nuclear interaction. In 
Figure [23] we give the case for the longest range. 
Here, large numbers of direct electron pair produc- 
tion with rather small energy loss play an impor- 
tant role, as shown similarly in Figure \T§\ Here, 
80.2% of the total energy is lost by 13722 direct elec- 
tron pair productions, 7.53% by 71 bremsstrahlungs 
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Figure 20: The energy losses with the shortest range for 10 eV muons. The figure is a magnification of Figure [2T1 
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Figure 21: The energy losses with the shortest range for 10 15 eV muons. 
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Figure 22: The energy losses with the average-like range for 10 15 eV : 
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Figure 23: The energy losses with the longest range for 10 eV muons. (For interpretation of the references to color in 
of these Figures 21 to 24, the reader is referred to the web version of this article.) 



40 

legend 



and 11.1% by 105 photo nuclear interactions. In 
Figure [22j we give the case with the average like 
range. Here, 23.5% of the total energy is lost by 
5489 direct electron pair productions, 23.5% by 49 
bremsstrahlungs and 0.93% by 37 photo nuclear in- 



teractions, while in the real averages, 47.4% of the 
total energy is lost by 6800 direct electron pair pro- 
ductions, 35.3% by 48.1 bremsstrahlung and 16.7% 
by 55.0 photo nuclear interactions. 

In Figures l24l to l27l we show the similar relations 
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Figure 24: The energy losses with the shortest range for 10 18 eV muons. The figure is a magnification of Figure [25] 
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Figure 25: The energy losses with the shortest range for 10 eV muons. 
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Figure 26: The energy losses with the average-like range for 10 °eV muons. 
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Figure 27: The energy losses with the longest range for 10 18 eV muons. (For interpretation of the references to color in '. 
of these Figures 25 to 28, the reader is referred to the web version of this article.) 



for 10 18 cV muons as shown in Figures [20] to [23] 
The case with shortest range in Figure [24] (Figure 
[231) has a strong contrast to that with the longest 
range. The manner of the energy loss in Figure [2"4l 
is drastic with two big catastrophic energy losses 



due to bremsstrahlung(~0.8 km, ~4.7 km), while 
that in Figure [27] is moderate with no catastrophic 
energy loss. The shortest range , ~7.7 kilometers, is 
far shorter compared with the longest range, ~57.8 
kilometers. It is seen from Figure [M] and the Ta- 
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Table 2: The details of the characteristics on the muons with the shortest range, the average-like range ,the longest range and 
the average range. 



E = 10 12 eV 


Range 
[km] 


Energy loss 
by brems 


Number of 
interaction 


Energy loss 
by direct pair 


Number of 
interaction 


Energy loss 
by nuclear 


Number of 
interaction 


< Average> 


2.43 


l.lOxlO 11 


4.74 


1. 57x10" 


243 


4.54xl0 lu 


3.44 


Shortest 


4.65xltr 2 


2.15X10 11 


6 


1.52X10 11 


208 


3.72xl0 IU 


3 


Average-like 


2.43 


8.97xl0 lu 


4 


1.34X10 11 


221 


8.86xl0 IU 


3 


Longest 


3.14 


3.80 xlO M 


5 


1.04X10 11 


299 


1.19xlO y 


3 


E = W 15 eV 




<Average> 


1.78X10 1 


3.53xl0 14 


48.1 


4.74xl0 14 


6.80xl0 :l 


1.67xl0 14 


55 


Shortest 


9.44X1C)- 1 


8.66xl0 14 


2 


1.34xl0 14 


367 


5.90xlO IU 


1 


Average-like 


1.78X10 1 


7.50xl0 14 


49 


2.35 xlO 14 


5489 


9.31xl0 12 


37 


Longest 


3.50X10 1 


7.53xlO la 


71 


8.02 xlO 14 


13722 


l.llxlO 14 


105 


E = 10 ls eV 




<Average> 


3.28X10 1 


3.37x10" 


108 


4.39x10" 


2.57xl0 4 


2.25x10" 


172 


Shortest 


7.72x10° 


8.75x10" 


28 


1.19x10" 


5760 


6.23xlO Ib 


40 


Average-like 


3.28X10 1 


1.68x10" 


118 


5. 58x10" 


29321 


2.74x10" 


196 


Longest 


5.78X10 1 


5.71xl0 1,J 


162 


5.77x10" 


46542 


3.66x10" 


277 



Table 3: The Ratios of energies transferred from 
bremsstrahlung, direct electron pair production and photo 
nuclear interaction to the total energy loss. 



E = 10 12 eV 


Brems 


Direct Pair 


Nuclear 


<Avcragc> 


3.37X10- 1 


5.26X10- 1 


1.37X10- 1 


Shortest 


5.32xT0 _1 


3.76X10- 1 


9.20xl0~ 2 


Average- like 


2.87X10- 1 


4.30X10" 1 


2.83X10- 1 


Longest 


3.50xl0- 2 


9.54X10- 1 


l.lOxlO^ 2 


E a = 10 lb eV 




<Avcragc> 


3.40X10- 1 


4.98X10- 1 


1.62X10- 1 


Shortest 


8.66X10- 1 


1.34X10- 1 


5.90xl0- 5 


Average- like 


7.54X10" 1 


2.37xi0 _1 


9.36xl0~ :i 


Longest 


7.62xl0~ 2 


S.llxlO- 1 


1.12X10- 1 


E„ = W ls eV 




<Avcragc> 


3.24X10" 1 


4.59X10" 1 


2.17xl0 _1 


Shortest 


8.75X10- 1 


1. 19x10^ 


6.23xl0~ ;i 


Average- like 


1.68X10- 1 


5.58X10- 1 


2.74X10- 1 


Longest 


5.71xl0~ 2 


5.77X10- 1 


3.66X10- 1 



bles that bremsstrahlung plays a decisive role as the 
cause of catastrophic energy loss. 87.5% of the to- 
tal energy is lost by 28 bremsstrahlungs, 11.9% by 
5760 direct electron pair productions and 0.623% 
by 40 photo nuclear interactions. In Figure 1271 we 
give the case with the longest range. Here, 57.7% 
of the total energy is lost by 46542 direct electron 
pair productions, 36.6% by 277 photo nuclear inter- 
actions and only 5.71% by 162 bremsstrahlungs in 
the complete absence of catastrophic energy losses. 
In Figure [551 we g ive the case average-like range. 
Here, 55.8% of the total energy is lost by direct elec- 



tron pair production, 16.8 % by 118 bremsstrahlung 
and 27.4% by 196 photo nuclear interactions, while, 
in the real averages, 43.9% of the total energy is lost 
by 2.57x 10 4 direct electron pair productions, 33.7% 
by 108 bremsstrahlung and 22.5% by 172 photo nu- 
clear interactions. Thus, it is can be concluded 
that the diversity among muon propagation with 
the same primary energy should be noticed. 

In Table [3j we give the proportion in the energy 
losses among bremsstrahlung, direct electron pair 
production and photo nuclear interaction for differ- 
ent ranges. It is easily understood from the table 
that in the cases with the shortest range, we obtain 
the proportions occupied by the bremsstrahlung 
5.32 x Mr 1 , 8.66 x 10" 1 and 8.75 x 10" 1 for 10 12 
eV, 10 15 eV, and 10 18 eV, respectively, while in 
the cases with the longest ranges, we obtain the 
proportions occupied by direct electron pair pro- 
duction, 9.54 x 10" 1 , 8.11 x 10" 1 , 5.77 x 10" 1 , re- 
spectively. Namely, catastrophic energy losses by 
bremsstrahlung make decrease their ranges remark- 
ably, while the direct electron pair production has 
high probability to lose smaller energy and this ef- 
fect makes it possible their ranges to be extended. 
Also, it is seen from the information around [aver- 
age/, that, on the average, a half of the total en- 
ergy loss is occupied by the direct electron the en- 
ergy loss and the ratio of the energy loss due to 
bremsstrahlung which is the origin of the strong 
fluctuation should be considered carefully. 
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3. Conclusion 

Many authors have examined the survival prob- 
ability on high energy muons by the Monte Carlo 
simulation which we call the v cut procedure, expect- 
ing to utilize their results for the analysis of the 
muon neutrino events in the KM3 detector in fu- 
ture. We have developed the time sequential pro- 
cedure in the Monte Carlo simulation without in- 
troducing v c ut) by which we give the same results 
with those by the v cu t procedure as far as the muon 
energy behavior is concerned. 

Using our method, we have examined the diver- 
sity on the muons' ranges (or, energy losses) in de- 
tail. Also, we have clarified that range energy distri- 
bution is well presented by the normal distribution 
over 10 14 eV to 10 18 eV and we have examined the 
degree on the fluctuation among three stochastic 
processes, by using the hypothetical range distri- 
butions due to the stochastic characters. 

Since we are only interested in behaviors of the 
high energy muon, there is no essential differences 
between the results obtained by the v cu t proce- 
dure and the time sequential procedure. However, 
the difference in the results obtained from both 
the time sequential procedure and the v cu t proce- 
dure may appear in Cherenkov light yields. As the 
Chercnkov light yields are closely related to the en- 
ergy estimation of the muon concerned, the differ- 
ence should be carefully examined. In a subsequent 
paper (Part2)[27]], we will discuss Cherenkov light 
yields due to high energy muon in greater detail 
with regard to muon straggling. 
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Appendix A. The Monte Carlo procedure 
for the time sequential proce- 
dure 

If the differential cross section for 
bremsstrahlung [28j . direct pair production(29| 
and photo nuclear interaction [30| are denoted by 
a b (E, Ej) cLE 7 , a d (E, E e ) dE e and a„ (E, E n ) dE n , 
respectively. Then, the mean free paths of the 
muons due to bremsstrahlung, direct electron pair 




Energy(GeV) 



Figure A.l: Mean free paths for bremsstrahlung (Curve b), 
direct electron pair production (Curve d), photo nuclear in- 
teraction (Curve n). 



production and photo nuclear interaction are given 
respectively as 



X d (E) 



A„ (E) 



a $ e iZZ ab ( E > E i) dE ~f 
i 

7T Ie"" 1 ^ °d (E,E e )dE e 



A JE n 



o~ n (E, E n ) dE v 



(A.l) 
(A.2) 
(A.3) 



Also, the resultant mean free path for these radia- 
tive processes are given as 

1 111 



XtotadE) X b (E) X d (E) X n (E) 



(A.4) 



The integrations for (A.l) to (A.3) are performed 
over kincmatically allowable ranges. In (A.l), 
E-y^nin is taken to be satisfied with such a con- 
dition that E lyTn i n / 'Eq is sufficiently smaller than 
Emin/Eo, where Eq and E m i n denote the primary 
energy of the muon and the minimum energy of 
the muon for observation. E m i n is taken as 1 GeV 
throughout present paper. In Figure A.l, we give 
the mean free paths for bremsstrahlung, direct elec- 
tron pair production and photo nuclear interaction 
are given as the function of the primary energy. 

The most important procedures in our Monte 
Carlo method are only two. The first one is where 
the interaction occurs. For bremsstrahlung, the tra- 
versed distance for the interaction is determined 
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Figure A. 2: Flow chart of muon propagation. 

with the use of £, a uniform random number be- 
tween (0,1), as follows. 

At = -A& (E) log£, (A.5) 
Similarly, for direct electron pair production, 

At = -X d (E) logt, (A.6) 
Similarly, for the photo nuclear interaction, 

At = -X n (E)log^ (A.7) 

In our Monte Carlo simulation, the energy losses 
due to each stochastic process are sampled by 
the following equations with the use of £, the 
uniform random number, between (0, 1) . For 
bremsstrahlung 



From |0]), 

Ey = F b (E, £) , 
For direct electron pair production 



Jp a d (E, E e ) dE e 



JeZZ Ud E ^ dEe 

From (TA~T0l) . 

E e = F d (E,0, 
For photo nuclear interaction 

/e„" a n(E,E n )dE n 

J^n m ax ^ ^ fifi^ 

From (TA~T2|) . 

E n = F n (E, £) , 



(A.9) 



(A.10) 



(A.11) 



(A.12) 



(A.13) 



Exactly speaking, the samplings in equations (A. 8), 
(A. 10) and (A.12) are carried out in v, the ra- 
tios of the sampled energies to their primary en- 
ergies, not the sampled energies themselves. For 
the discrimination among stochastic processes in 
our Monte Carlo simulation let us introduce the 
following equations. 



ta(E) 



1/A& (E) 



l/X b (E) + 1/X d (E) + I/An (E) 

1/X b (E) + l/X d (E) 
l/X b {E) + l/X d (E) + 1/X n (E) 



j E ' 1 <J b (E , Ey) (lEry 

/ £ 1m " (J b (E, Ey) dEry 



(A. 



In Figure IA.21 we give a flow chart for our Monte 
Carlo simulation. The validity of our Monte Carlo 
procedure by the time sequential procedure should 
be carefully examined in the two methods which 
are independent of each other. Namely, the first 
one is the comparison of our procedure with the 
analytical theory which is explained in the present 
Appendix and the other is the comparison of our 
procedure (the time sequential procedure) with the 
different kind of Monte Carlo procedure (the v cu t 
procedure) which is mentioned in the text. Par- 
ticularly, it is best that the results obtained by a 
Monte Carlo procedure are checked by the proce- 
dures which are methodologically independent of 
the Monte Carlo procedure, reaching the same re- 
sults. In Figure A. 3, the average energies of the 
muons are given as a function of the depths un- 
der the preposition of muon energy spectrum at sea 



13 




Figure A. 3: The average energies of the muons. The lines 
denote our, while symbols denote Misaki and Nishimura. 



level with indexes 2, 3 and 4, obtained by the time 
sequential procedure [2~i| p6| and they are compared 
with results obtained by the analytical theory based 
on the Nishimura-Kamata formalism in the cascade 
shower theory [III ] and the agreement between them 
is very good, which surely guarantees the validity 
of our time sequential procedure. 
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